function theta=computeTheta(identityC, multiColumnPowersAB)
theta = identityC * multiColumnPowersAB;
end